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Abstract. We report recent results of strangeness and heavy flavor measurements 
from PHENIX. The topics are: Elliptic flow of strangeness and heavy flavor electron 
production comparing to the other hadrons, </> meson production, and an exotic particle 
search. 

(Some figures in this article are in color only in the electronic version) 



1. Introduction 

This talk will focus on strangeness and heavy flavor and report recent results from 
PHENIX. The topics are the following: (1) event anisotropy of charged n, K, p and p 
in y/SNN = 62.4 GeV Au-|-Au collisions comparing to 200 GeV, (2) event anisotropy 
of heavy flavor electron in 200 GeV Au-|-Au, (3) the invariant cross-section of meson 
production in 200 GeV d-f Au and Au-|-Au, (4) exotic particle search in 200 GeV p+p, 
d-|-Au and Au-|-Au, and (5) detector upgrade for future strangeness and heavy flavor 
measurement at PHENIX. A detailed report and discussion about heavy flavor results 
is presented in jT]. 

The PHENIX detector consists of two central arms and forward muon arms. The 
central arms have the capability of identifying and measuring charged hadrons (vr^, 

p, p, d, and d), electrons and positrons, photons, and anti-neutrons. By using 
the invariant mass spectrum, weak-decay particles and resonances can be identifled, 
for example, 7r°, cj), A (A), and E^. The forward muon arms identify /i^. A detailed 
description of the PHENIX detector can be found in Pj. 



2. Event anisotropy of hadrons 

Event anisotropy analysis is a powerful tool to study properties of the early stage in high- 
energy heavy ion collisions. The second Fourier coeflicient of the azimuthal momentum 
distribution, t>2, characterizes the event anisotropy. Recent results of charged vr, K, p, 

I For the full PHENIX collaboration author list and acknowledgments, see appendix 'Collaboration' 
of this volume. 
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Figure 1. (Color online) The left plot shows PHENIX preliminary results of charged 
TT, K, and p U2 as a function of pr from minimum bias data in ^smn = 62.4 GeV 
Au+Au collisions. The tt^'^, 7r°, p, and p V2 in 200GeV Au+Au collisions are 
shown in the right-hand side. The identified charged hadron V2 is from W and tt" data 
has been presented in The systematic errors related to particle identification and 
background are similar in both 62.4 GeV and 200 GeV. The difference comes from 
systematic uncertainties in the reaction plane determination due to different statistics 
and multiplicity in the BBC between 62.4 GeV and 200 GeV. 



and p V2 are shown in figure Q in 62.4 GeV and 200 GeV Au+Au collisions. Both 62.4 
GeV and 200 GeV results are from minimum bias trigger data (about 40M and 30M 
events, respectively). The charged hadrons are identified by Time-Of-Flight (TOF) 
information with momentum. The neutral pion is reconstructed from photon pairs 
measured by the Electro-Magnetic Calorimeter (EMCal). The V2 is computed with 
respect to the reaction plane which is defined by hits in the Beam-Beam counters (BBC) 
located in the forward/backward regions (north and south) at |?7| = 3.1-3.9. Since the 
V2 measured is smeared due to reaction plane resolution, it is necessary to correct for 
their effect. The correction factor is calculated from the two reaction planes defined by 
data from the north and south BBC separately, using a method described in 0. 

The V2 in i/Jntv = 62.4 GeV Au+Au collisions shows the same order of magnitude 
and mass dependence as in 200 GeV Au+Au. As we reported in H], there is a 
crossover point among hadrons around pt = 2GeV/c in 200 GeV Au+Au. The proton 
V2 crosses the pion and kaon around pr = 2 GeV/c in 62.4 GeV Au+Au, the same as 
in 200 GeV. Since the lower p^ (< 2 GeV/c) region in 62.4 GeV Au+Au shows similar 
mass dependence to 200 GeV Au+Au, it might be well described by hydrodynamical 
models. In the pr>2 GeV/c, the meson V2 remains different from that for baryons. 
This suggests a scaling rule from a quark coalescence model jB], which relates to the 
different number of quarks in mesons and baryons. Additionally, it might indicate that 
the hadron f 2 is generated during the early partonic stage before the hadronization. 

FigureElis a simple check of the coalescence picture in f 2 for normalizing both V2 and 
Pt by the number of constituent quarks. The data in 62.4 (200) GeV Au+Au collisions 
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Figure 2. (Color online) V2 as a function of px scaled via the coalescence prescription. 
The solid (open) markers show the data from ^/snn = 62.4 (200) GeV Au+Au. 



are shown by solid (open) markers in the figure. We make the following observations: 
(1) The 62.4 GeV data seem to be described by the coalescence picture, (2) the 62.4 GeV 
data are slightly smaller than 200 GeV but consistent within the systematic uncertainty, 
and (3) observables (1) and (2) suggest a saturation of azimuthal anisotropy in Au+Au 
collisions in the range from 62.4 to 200 GeV. More statistics in 62.4 GeV Au+Au is 
needed to conclude whether the particle type dependence of as a function pt and 
centrality are the same as in 200 GeV Au+Au collisions. 

3. Event anisotropy of "non-photonic" electrons 

PHENIX has measured inclusive electrons (and positrons). The data include two 
components, which are 'photonic' and 'non-photonic' electrons. The 'photonic' electrons 
mainly come from Dalitz decay (of vr", 77, cu, 0, and etc.), di-electron decays (of p, cu, 0, 
and etc.), photon conversions, and kaon (^ vrez/) decays. The 'non-photonic' electrons 
are estimated by subtracting 'photonic' electrons from inclusive electrons. (A detailed 
discussion of the electron measurement is found in reference 0.) Although the non- 
photonic electrons are decay products from both charmed and bottom mesons, event 
anisotropy for the pt range from the Run-2 data set are dominantly charmed meson 
decays. 

The 'non-photonic' electron V2 in 200 GeV Au+Au collisions from Run-2 is shown 
in figure El with statistical error (vertical bar) and systematic errors (shaded box). The 
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Figure 3. (Color online) "Non-photonic" electron V2 as a function of pr in ^/s^n = 
200 GeV minimum bias Au+Au collisions from Run-2. The vertical bar on each point 
shows the statistical error and systematic errors are shown by the colored box. The 
horizontal bar shows the RMS of dN/dtpx in each bin. Hash marks at the top indicate 
bin boundaries for the pr bins and the last bin is an average of pt = 1-75 to 4.0 
GeV/c. The solid and dashed lines are model calculations shown in [7]. The solid line 
shows the prediction for D meson decays which are completely thermalized, including 
transverse expansion for the charm quark. The dashed line is for D decays which have 
no rescattering, corresponding to perturbative QCD spectra for the c-quark. 

two lines in figure 01 address two scenarios from a model calculation in The solid line 
shows one scenario in which D mesons are made from completely thermalized charm 
and light quarks; the other (dashed line) is that there is no interaction of the c-quarks, 
so that the flow contribution in D's is only from the light quarks. Within the available 
statistics, both models are consistent with data. It will be very interesting to pursue 
this comparison of charm flow with that of mesons containing light quarks using the 
higher statistics Run-4 data set. Additionally, the systematic error of the 'non-photonic' 
electron V2 is expected to be relatively smaller for the Run-4 data analysis, because 
uncertainty of the cross-section and V2 for 'photonic' electrons will be reduced. 

4. meson production 

The meson invariant cross sections in ^snn = 200 GeV d-l-Au and Au+Au collisions 
are shown as a function of mT-m^ in figure 01 where is the transverse mass 
{=\Jpr^~-\-^rn^) and is the mass from PDG values. In d+Au collisions, 31M 
single-electron triggered events were used for the e~^e~ analysis and 62M minimum bias 
events were used for the K^K~ . 20M minimum bias triggered events are used in the 
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Figure 4. (Color online) The left-hand plot shows the invariant cross-section for (j) 
production as a function of (m-r - mass of </)) for the indicated centrality bins in ^/snn 
= 200 GeV Au-I-Au colhsions (from [S]). The data points are for the K+K~ 
channel only. The lines are my exponential fits (see text in detail). The results from 
y/SNN = 200 GeV d+Au collisions are shown in the right-hand. Statistical error bars 
are shown. 




analysis of Au-|-Au collisions. These numbers include a vertex cut (l-Zvertoxl < 30 cm). 
The spectrum is fit with the exponential function: 

2'Kmj' dmrpdy 2'kT(T -|- m^) 

where dN/dy and the inverse slope parameter T are fit parameters. The fit results are 
shown in figure El and table ^ The left-hand plot of figure El show the inverse slope 
parameter as a function of average number of participants {{Npart)) in 200 GeV d-l-Au 
and Au-|-Au collisions. The dN/dy normalized by (Npart) are shown in the right-hand 
of figure The Au-|-Au results are from the (p —>■ K^K~ channel only, while results 
from — i> K^K~ and e^e~ channels are presented in d-l-Au collisions. In the 
d-|-Au collisions, the invariant cross sections from the two channels show agreement 
around tut — = 1 GeV/c. Both the inverse slope parameter and the yield for the 
two channels are consistent, as can be seen in figure 

Because of large statistical and systematic errors in the electron channel, (Npart) 
dependence of phi meson will be discussed with only the kaon channel data in the 
following. The inverse slope parameter as a function of (Npart) seems to be fiat over 
all (Npart) bins within errors in d-|-Au and Au-|-Au collisions. On the other hand, 
[dN/dy)/ (Npart) increases by about a factor of two from d-l-Au to Au-|-Au data. 
Within the statistical and systematic error bars, the yield normalized to the number of 
participants is independent of centrality in Au-fAu collisions. 

The ratio Rcp (central to peripheral ratio scaled by number of collisions) of protons 
{{p+p)/'2) is different from neutral pions as a function of pt in ^/SNN = 200 GeV Au-|-Au 
collisions. As shown in figure IHl proton Rcp increases up to = 2 GeV/c and appear 
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Figure 5. (Color online) The inverse slope parameters (left) and the yield normalized 
by the average number of participants {Npart) (right) of (j> in ^snn = 200 GeV d+Au 
and Au+Au collisions as a function of (Npart) are shown. Those values are obtained by 
a fit using equation^ The fits are done separately for the electron and kaon channels. 
The vertical bars show the statistical error only. The systematic errors are shown as 
shaded boxes. 



Table 1. The inverse slope parameters and the yield of 4> meson in ^/snn = 200 GeV 
d+Au and Au+Au collisions. Values are obtained using equation^ The fits are done 
separately for the electron and kaon channels. (Npart) is shown with only statistical 
errors. Errors of the inverse slope and the yield are statistical (first) and systematic 
(second) . 



Collision Decay Inverse slope 

system (Npart) mode [MeV] dN/dy 



d+Au (Min. Bias) 


9.2+0.8 




326 + 94+173 


0.056 + 0.015 + 0.028 






K+K- 


414 + 31 + 23 


0.0468 + 0.0092+°;°°^^ 


Au+Au (40-92%) 


32.0+2.9 


K+K- 


359+15 + 16 


0.32 + 0.03 + 0.05 


Au+Au (10-40%) 


171.8+4.8 


K+K- 


360+13 + 23 


2.22 + 0.18 + 0.35 


Au+Au (top 10%) 


325.2+3.3 


K+K- 


376+24 + 20 


3.94 + 0.60 + 0.62 



to saturate around 1, up to px = 4.5 GeV/c. The 7r° Rcp is about 0.4 over pt= 1 to 
7 GeV/c. The (p Rcp is about 0.6 over px = 0.8 to 3.0 GeV/c and closer to the pion 
than proton. It suggests that there is a scahng with the number of valence quarks in 
Rcp- On the other hand, a detailed discussion on such scaling should be based on more 
hadrons, for example, Kg, A, S, H, and Q. 

5. Exotic particle search 

One of the most interesting and important topics in hadron physics is the possible 
existence of a five-quark state. The LEPS collaboration has reported evidence for a 
narrow 5=+! baryon resonance state ^U]- The reported resonance is a candidate 
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Figure 6. (Color online) Rcp (central to peripheral ratio scaled by number of 
collisions Ncou for 0, {p + p)/2, and 7r° are shown as a function of in ^/snn = 200 
GeV Au+Au collisions. The shaded solid bar around Rcp=l is 12% systematic error of 
the ratio. The vertical dotted bar on the right presents the error on N^oii^^°/'^coii^^'^°- 
The dotted horizontal line on the cj) data points is a straight-line fit to the (/) data. 



for the state of Q^{uudds) which was predicted with mass 1.540 GeV. has two 

decay modes, K^+n and K^+p. Additionally, the NA49 experiment at CERN-SPS 

has reported additional potential pentaquark states that are called E'^^ {ssddu) and S3 

22 

(ssudd) p.!] . 

Since the anti-baryon to baryon ratio is larger than 0.75 at mid-rapidity in ^/snn 
= 200GeV p+p and Au+Au collisions jT^ IT^ IT^ [T3] , we expect a similar or larger anti- 
pentaquark to pentaquark ratio than 0.75. We search for the pentaquark by using an 
anti- neutron, that is, Q~^K~+n. The advantages are the following. The anti- neutron 
will be identified by the Electro-Magnetic Calorimeter (EMCal), and the K~+n channel 
(two-body decay) has a larger acceptance than Kg+p{p) (with three final state, i.e., 
TT~^+7r~+p (p)). 

The anti-neutron leaves a large energy deposit and a different energy shower shape 
than for photons and electrons. Therefore, we identified anti-neutrons requiring three 
conditions: (1) An energy deposit in EMC is larger than 1 GeV, (2) the energy deposit 
is not due to a charged particle, (3) the shower shape is different from photons and 
electrons. The neutron momentum is reconstructed by the TOP in the EMcal and path 
length. The momentum range of the neutrons is 0.5 to 3.0 GeV/c in px, with the above 
conditions. The K~ is identified by a cut of momentum and EMCal Time-of-Plight 
(within 3 0" of TOP resolution, with 1 a about 500 ps). 

The resolution of momentum calculated by EMCal TOP is studied with an anti- 
proton sample. We compared momentum from two calculations, a normal tracking 
method for charged particles and the other from TOP. The results shows that the 
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resolution is less than 4%, 10%, 20% in ]9r<1.0, 1.4, 3.0 GeV/c, respectively. Another 
study was done for invariant mass for (tt^ + n) reconstruction. The peak was 
confirmed at the expected mass within 10 MeV for pt<3 GeV/c, and the results 
were consistent with the momentum resolution study. 

The data set of minimum bias trigger that we used consists of 35M, 91M and 36M 
event in ^Jsnn = 200GeV p+p, d+Au and Au+Au collisions, respectively. The invariant 
mass distribution of K~ and n pair is calculated in each event. The distribution has 
a combinatorial background and we estimate it by event mixing, i.e., by pairing tracks 
of similar vertex and centrality from different events. The results show no significant 
pentaquark peak around expected mass (~ 1.540 GeV/c^) in p+p, d+Au and Au+Au 
collisions. 

6. Detector upgrade 

There are several projects to enhance the particle identification capability for the higher 
momentum region, and to improve the momentum resolution in PHENIX. Aerogel 
Cherenkov Counters have been installed to cover |r7|<0.35 and A</)=7r/8 in 2003 and 
2004. Proton separation from tt/K will be done up to p=7 GeV/c (The current 
maximum momentum is 4 GeV/c by TOF detector). The following four detectors 
are planned. A silicon vertex detector and TPC are planned to enhance momentum 
resolution and to measure in-fiight-decay particles. They will be located close to the 
collision vertex and will reduce background in the high pt region from in-fiight-decay. 
A hadron blind detector will surround the TPC and has the capability of electron ID by 
Cherenkov radiation detected by a photocathode layer. To measure jets in the forward 
region, a nosecone calorimeter is being designed and will cover ?7=0.7-2.6. 
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